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Abstract 

Aims. We study the structure of the medium surrounding sites of high-mass star formation to determine the interrelation between the H ii regions 
and the environment from which they were formed. The density distribution of the surroundings is key in determining how the radiation of the 
newly formed stars interacts with the surrounds in a way that allows it to be used as a star formation tracer. 

Methods. We present new Herscliel/SPIRE 250, 350 and 500 yum data of LHA 120-N44 and LHA 120-N63 in the LMC. We construct average 
spectral energy distributions (SEDs) for annuli centered on the IR bright part of the star formation sites. The annuli cover ~ 10 — 100 pc. We use a 
phenomenological dust model to fit these SEDs to derive the dust column densities, characterise the incident radiation field and the abundance of 
polycyclic aromatic hydrocarbon molecules. We see a factor 5 decrease in the radiation field energy density as a function of radial distance around 
N63. N44 does not show a systematic trend. We construct a simple geometrical model to derive the 3-D density profile of the surroundings of these 
two regions. 

Results. Herschel/SPIRE data have proven very efficient in deriving the dust mass distribution. We find that the radiation field in the two sources 
behaves very differently. N63 is more or less spherically symmetric and the average radiation field drops with distance. N44 shows no system- 
atic decrease of the radiation intensity which is probably due to the inhomogeneity of the surrounding molecular material and to the complex 
distribution of several star forming clusters in the region. 
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1. Introduction 

High-mass star formation (SF) sites (hereafter HMSFSs) are the 
beacons by which we probe a large part of the physics of exter- 
nal galaxies. They generally represent the most important tracers 
of the properties of their host galaxies in terms of star formation 
rate (SFR) and general activity. The most frequent tracers of the 
star-formation activity generally use the fact that the abundant 
UV light coming from the hot, young stars is absorbed in the 
vicinity and reradiated in the form of line or continuum emis- 
sion. This is true, for example, for Ha (e.g. Kennicutt,1998^, the 
aromatic emission bands in the mid-IR (e.g. Calzetti et al.|2007"]l 
or the IR continuum due to solid-state materials (dust, e.g. fSoifer 
let al.||1986|l. These tracers work relatively well and are used to 
characterise nearby star forming regions and star forming galax- 
ies out to large redshifts, although interesting discrepancies have 
been noted for dwarf galaxies at low SFR (see Lee et al.|2009 1. 
One of the main assumptions that enters into the quantitative 
interpretation of these data is the geometry of the material sur- 
rounding the newly formed stars, in particular, where the UV 
light is being reprocessed. For example, if the UV photons es- 
cape from the ionised medium this may boost the aromatic fea- 
ture strengths and strongly influence the lines originating from 
the surrounding photo-dissociation regions (PDRs). There are 
indications that this geometry in external galaxies may quali- 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



tatively and quantitatively differ from that observed in the MiUcy 
Way (MW) (e.g. |Madden et al.|2006[|Galametz et al.|2009[ ). One 
simple effect may be that at different metallicities the surround- 
ing medium is more or less opaque and therefore the UV pho- 
tons have a different mean free path. More complex scenarios are 
also discussed in the literature. For example, dumpiness of the 
molecular cloud may lead to small molecular cores surrounded 
by large PDRs. Herschel with its unprecedented wavelength cov- 
erage and angular resolution at submillimeter (submm) wave- 
lengths provides a unique opportunity to probe the cold inter- 
stellar medium (ISM) and sample the effects of the environment 
on the resulting SF tracers. In particular, it is well suited to trace 
the distribution of matter around HMSFSs and to map the way 
the UV radiation permeates and heats the surroundings. Here 
we present a study of two HMSFSs in the Large Magellanic 
Cloud (LMC) ba sed on data taken in the HERITAGE program 
(PI. Meixner, see |Meixneret al. 2010). 

2. Data treatment and modeling 

We construct spectral energy distributions (SEDs) of the envi- 
ronment of two distinct HMSFSs (LHA 120-N 44 and LHA 120- 
N 63 (Henize 1956 ), hereafter N44 and N63, respectively) in the 
strip of the LMC that was mapped during the science demonstra- 
tion phase. The two regions were chosen because they are the 
brightest and most isolated regions observed. N44 is the bright- 
est H II complex in the observed strip. Massive star formation 
has been active in this region; it contains three OB associations 
LH47, 48, and 49 with ages > 10 Myr at the central super-bubble 
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Figure 1. False colour images of the surround- 
ings of N44 and N63. IRAC 4 (8.0 fim; blue) 
traces the PAH emission, i.e. dominated by 
the PDRs that are the result of the UV radia- 
tion impinging on the surfaces of the molec- 
ular clouds. MIPS 24 fira (green) traces the 
hottest dust and is generally a good tracer of 
the ionised medium around SF sites. SPIRE 
250 fim (red) maps to first order the dust col- 
umn density and is less sensitive to the ac- 
tual illumination conditions. The circles desig- 
nate the annuli for which we construct SEDs. 
The coordinates of the centers are 05:22:03.2 - 
67:57:55 and 05:35:41.3 -66:01:45 [J2000] for 
N44 and N63, respectively. The yellow circles 
indicate the positions of optically detected ion- 
ising stars jOey & Massey|1995[[Oey|1996[ > 



and < 5 Myr at the super-bubble rims and surrounding dense 
H II regions, as well as a large number of massive young stel- 
lar objects (Y SOs) with ages < 1 Myr ( |Lucke & Hodge||1970[ 
Oey & Massey 1995 , Chen et al. 2009). N44 is also the brightest 
source at al l of the SPIRE wavelengths (250, 350, and 500 //m) 
in the strip ( Meixner et al.,2010^ . By contrast, N63 is a simple, 
roughly round H ii region. It contains one OB association LH83 



with an age < 5 Myr and a number of massive YSOs ( Bica et al. 
T996l|CauIetetaI.|2008| ). At SPIRE wavelengths, N63 is in a rel- 



atively isolated environment. Thus, we use N44 as an example of 
a prominent H ii regions seen in more distant galaxies and N63 
as a comparison whose simple structure makes it more straight- 
forward to relate dust properties with physical conditions of the 
ISM. 

The data we use are: 2MASS J, H, K, fSk rutskie et al. 2006), 
Spitzer I RACl, 2, 3 4 ([Fazio et al.||20()4|, MIPS 24, 70 and 
160 fim ( |Rieke et al.||2004| l a nd Herschel |pabratt et aI.|[20T0ll 



SPIRE 250, 350 and 500 yum ( jGrifRn et al. 2010|l. See [Meixner 



et al.[ ( [2010^ for a description of the SPIRE data treatment. We 
do not use the PACS information since the data we have until 
now does not allow us to extract extended source fluxes with 
sufficient accuracy. We have extracted maps of 40' x40' centered 
on each H ii region. We have convolved these data to a spatial 
resolution of 38" set by the MIPS 160 jum/SPIRE 500 /urn data. 
The 2MASS data have been convolved with the beam of MIPS 
1 60 fim, the IRAC 1 -4, MIPS 24, 70 /vm have been convolved us- 
ing custom-made kernels (Gordon et al. 2008 ). The SPIRE data 
have all been convolved to the SPIRE 500 fim resolution assum- 
ing Gaussian beam profiles withFWHMof 18.1, 25.2 and 36.9" 
for SPIRE 250, 350 and 500, respectively. We have resampled 
the convolved images to the pixel scheme of the MIPS 160 fim 
image using the IDL/astrolib routine hastrom. 

The center of the HMSFS is determined, by fitting a 2D 
Gaussian profile, as the brightest source in the total IR (TIR) 
image. The TIR image is obtained by simple integration from 
8 to 500 /vm. The reasoning behind using this definition of the 
center is that this location (on a size scale of tens of parsecs ) 
probably hosts the most active site of embedded star formation. 
We extract flux densities (in Jy) for annuli around the given cen- 
ter (see Fig. [Til. We use the following radii for the sizes of the 
annuli: 35, 47, 78, 125, 200, 330 and 530", corresponding to lin- 
ear sizes of 8, 12, 19, 30, 48, 80 and 128 pc assuming a distance 
to the LMC of 50 kpc (e.g. [Schaefer|[2U08) . The annuli were 
chosen to be larger than the apparent size of the HMSFS so as to 
be able to also study the regime in which the emission becomes 
dominated by the general LMC. Examples of the extracted SEDs 
for each region in different annuli are shown in Fig.|2^,b. There 



is a clear and systematic trend for the far-IR to peak at longer 
wavelength for the outer annuli. In the case of N44 this trend is 
to a large part offset by the broadness of the far-IR peak. This is 
an indication that the emission arises from a broad temperature 
distribution. We also note that we detect polycyclic aromatic hy- 
drocarbon (PAH) emission (at 8 //m) in all SEDs. Since we are, 
to first order, interested in the "shape" of the IR SED as a func- 
tion of distance, we create relative flux density maps by dividing 
the maps at each wavelength by the TIR map in order to measure 
the scatter. The scatter inside each annulus of the relative inten- 
sity map was used to estimate the uncertainty (Icr) on the flux 
densities. Note, that this method for determining the variance 
causes a large scatter in the near-IR pixels which are dominated 
by stellar light. This is the reason for the large error-bars at these 
wavelengths in Fig.[2^,b. 

We use a phenomenological dust model, which is described 
( |2009l l. 



in detail in Galametz et al. 



to interpret the observed 
SEDs. This model consists of a quantity of dust, with a realistic 
composition (PAHs, silicate and graphite) and grain-size distri- 
bution (complex molecules, very small (VSG) and big grains) 
illuminated by a radiation field with a range of intensities. The 
main parameters that we aim to constrain are: 1) total dust mass, 
2) the range of illumination intensities, and 3) the relative com- 
position, in particular the mass fraction of PAHs. The radiation 
field distribution is represented by a power-law function that de- 
scribes which fraction of the matter is exposed to what radiation 
field Pale & H elou 2002): 
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where U is the intensity of the incident radiation field the dust is 
subjected to (U-\ for the diffuse MW environment), [/,„!„ and 
L'^max the minimum and maximum value of U and <U> is the 
mass-weighted averaged radiation field intensity. 

The starting point for the modeling is to adopt dust proper- 
ties that fit the MW IR emission well ( [Zubko et aLl|20()4[ dis- 
tribution BARE-GR-S). The observed 24 yum fluxes in the dif- 
fuse ISM force us to use a dust size distribution that is different 
from the MW values. The best fits are obtained by reducing the 
mass fraction of VSGs to total dust by 50 per cent. We use the 
MW interstellar radiation field as the shape of the radiation field 
and do not vary this. The derived dust masses are robust against 
the choice of interstellar radiation field and the mass fraction of 
VSGs. However, the relative mass fractions of the various dust 
constituents may depend on these choices. 
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Figure!. Summary of main 
results. SEDs of the two SF 
regions we study (left). Tlie 
black symbols show the inte- 
grated photometry of three dif- 
ferent annuli, lines are the best 
fit models and the coloured 
triangles represent the syn- 
thetic photometry in the cor- 
responding filters. The derived 
dust column density (panel c), 
mean radiation field energy 
density (panel d), and the 
PAH mass fraction (panel e) 
as a function of annulus size 
are shown on the right. The 
error-bars on the parameters in 
Panels c,d,e have been derived 
by propagating the variance on 
the photometry using a Monte- 
Carlo method (see Sec.[2]l. 



In comparing the dust mass tracers with the gas tracers 
HMeixner et al.|2010HGordon et al.|2010]|Duval et al.|2010| i some 
issues have been raised about the applicability of this composi- 
tion to the LMC, in particular the graphitic component. We have 
verified that the derived mass and radiation field profiles (the 
shape as a function of annulus) are not sensitive to the choice of 
the carbon bearing grains. However, the absolute values of the 
derived parameters depend on this choice. Uncertainties on the 
derived parameters were estimated using a Monte-Carlo evalua- 
tion. The fitting routine was repeated 300 times with the obser- 
vational constraints varied randomly, according to their standard 
deviations and new best fit parameters are determined. The en- 
semble of best fit parameters is used to calculate the error-bars 
on each parameter (see Fig. [2]). 

3. Results 

Fig.l2]a,b show representative SEDs and the best fit models. The 
models fit the data very well over the entire wavelength range, 
for all annuli. Most parameters are well constrained with a distri- 
bution of best fit parameters which is roughly symmetric around 
the central value. One exception to this is the <U> in the two 
outer annuli of N44. We find a distribution around the mean 
which is heavily skewed to low values of <U>. We do not detect 
a systematic submm excess, i.e. the whole wavelength range up 
to 500 fj.m is well fit by the standard model. The model and the 
500 fxm surface brightness agree on average within 3% with a 
very small dispersion. The model is in accordance with the find- 
ings of Gordon et al. (2010 1, who show that the evidence in the 
LMC for a submm excess is weak and if present it is confined to 
the more diffuse and fainter environments. 

The SED fits allow us to study the derived properties as a 
function of annulus (projected distance). We show in Fig.l2t,d,e 
as a function of radius, the dust column density (Sdust), average 
radiation field intensity (<U>) and the mass fraction of PAHs 
normalised to the MW value. Both sources show a column den- 
sity profile that decreases outwards, although the column density 



in the central region of N63 is not much higher than in the outer 
annuli. N63 causes an increase of a factor 2-3 of column den- 
sity. The column density towards the central regions of N44 is 
~10 times higher than its environment. 

The width of the IR SED, which is well determined using 
the SPIRE data, requires a significant dust mass at lower tem- 
peratures. This is reflected in the low values of <U> in panel d. 
<U> exhibits only a small range for the entire sample of SEDs. 
<U> decreases as a function of distance for N63. N44 does not 
exhibit any significant trend and the data are consistent with a 
constant <U> over a distance scale from 10-100 pc. Note that 
<U> is weighted by dust mass and so it readily traces the radia- 
tion field as seen by the coldest dust along the line of sight. 

We calculate the second moment of the radiation field distri- 
bution (At/), again weighted by dust mass in order to quantify 
the range of radiation field that the matter is exposed to. All SEDs 
require a wide range of U. For N44 AC/ is roughly constant at the 
value of 30 meaning that the entire region is typified by a radia- 
tion field intensity ranging from 1 to 30. N63 shows an outward 
decreasing At/ from ~50 in the inner annulus to 2 in the outer. 
Thus the radiation field in N63 spans 5-55 (<\J>-<\J>+AU) on 
the inside and 1-3 in the outer annuli. This implies that the outer 
annuli in N63 closely resemble a diffuse environment. 

The fraction of mass contained in PAHs increases with in- 
creasing radius out to about 50 per cent of the MW value. We 
find a significant depletion of the PAHs towards the central re- 
gions over ~20 and 40 pc, for N44 and N63 respectively. Note 
that the smaller size of the depleted region is consistent with the 
fact that N63 is less prominent compared to its surroundings (see 
below). Interestingly, the radius of the depleted environment cor- 
responds well to the area occupied by the ionising stars in Fig.fT] 
The fraction of the ionised PAHs is not well constrained but we 
have checked that there is no systematic effect of the fitting pro- 
cedure that causes the observed trend in Fig. l2fe. 

We construct a simple geometrical model to try to constrain 
the size and density of the matter in the HMSFSs, i.e. to depro- 
ject the observed column densities. We represent the molecu- 
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Figure 3. Derived geometrical parameters of the two regions. We 
show the dust column density and the fitted profile of a thick 
shell, with a power-law density profile (blue line) and a slab 
(red). The dashed model corresponds to the case when we fix 
the exponent of the power-law to 2 (see text). 



lar cloud around the H ii region as a geometrically thick shell. 
To this we add a slab that represents the disk of the LMC, i.e 
the diff'use extended emission. We assume a constant dust col- 
umn density for the slab and a power-law profile for the shell 

(Pdust = Pdust,in('"/^in)"'')- 

The dust mass profile and the best fit decomposition into 
these two components and the corresponding parameters are 
shown in Fig. [3] We have propagated the uncertainties on the 
dust column density in each annulus (see Fig. l3]l. The uncer- 
tainties are significant. We find in case N63, for which the over- 
density is weak compared to the slab (i.e. the disk of the LMC), 
that the geometrical parameters are not very well constrained. 
In particular, it is hard to constrain the y parameter. The best fit 
model has y - 1.2 which seems flat. If we fix y at 2, the value 
expected for a cloud in hydrostatic equilibrium, we find that the 
distance where the slab component starts to dominate, i.e. the 
intersection of the blue and red lines in Fig.lS] reduces from ~80 
to -30 pc. 

4. Discussion & Conclusion 

Figs.|2]and[3]show the power of Herschel to determine the matter 
distribution around HMSFSs. We investigate the effect the new 
SPIRE constraints on the derived parameters and their uncertain- 
ties by also fitting the SEDs without using these SPIRE data. We 
find that the column density of dust is often very discrepant (by 
more than an order of magnitude) from the values derived using 
the SPIRE data. The derived Z^,,,, profile for N44 (Uke Fig. |f] 
top panel) is virtually constant at a value of 0. 1 Mq arcsec"^ due 
to the fact that the cooler dust, is not well traced by the MIPS 
data. For N63 the model over-predicts the derived dust masses 
for the annuli between 20 to 60 pc by a large factor (~5). In these 
cases the model tends to find very small values of C/niin- This 
artifact could perhaps be circumvented, when modelling those 
kind of regions, in cases where submm constraints are missing, 
by limiting the allowed range of J/min in Eq. (fTl). The difficulty 
of constraining the radiation field parameters without the SPIRE 
bands is also reflected in the derived uncertainties. In particular 
the value of f/min is ill-constrained which results in large, asym- 
metric error-bars on the derived column densities. 



N44 and N63 exhibit strikingly different behaviour in the 
radiation intensity profile (Fig. |2ll.) The lack of a systematic 
decrease of <U> around N44 indicates that we are observing 
dust with a wide range of temperatures along each line-of-sight. 
The inner annuli in N44 are affected by the superbubble to the 
NE of the OB association, where high values of <U> are ex- 
pected. The low values of <U> for such a luminous SF region 
may reflect dumpiness. The profile is clearly incompatible with 
a centrally illuminated optically thin irradiation profile. It is clear 
from Fig.[T|that the studied regions are not very spherically sym- 
metric (azimuthally smooth). In particular, N44 harbours several 
clusters and the peak of the X-ray emission is located in a cavity, 
~ 20-30 pc away from the TIR peak (see Fig.[T]i. Measuring the 
azimuthally averaged properties smears out some of the charac- 
teristics. This smearing could have been the cause for the lack of 
trend seen in the average U as seen by the dust (Fig. [2]). We have 
verified that this small range of <U> is not simply an artifact of 
this averaging or a wrong choice of center of the annuli by study- 
ing the parameters we derive pixel by pixel in the maps which 
makes no assumptions about the geometry. Indeed the highest 
<U> in N44 is found close to the center we chose. Except for 
the very center all other values, with their scatter, are within the 
range as depicted in Fig.l2] We conclude that the choice of center 
does not dominate the lack of trend of N44 in the average radi- 
ation field. Thus this lack of trend reflects the true broad range 
of irradiation conditions along all lines of sight in N44 , which 
is an indication of the inhomogeneity of the ISM around N44. 
A simple dust model shows a deficit in PAHs toward the centers 
of these two regions. We find no evidence for a submm excess. 
We have used the observed dust column densities surrounding 
N44 and N63 to derive a 3-D model for these regions for the 
first time. 
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